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CN Abstract 

The spectrum of tetraquark states with hidden charm is studied within an effective 
Qh QCD Coulomb gauge Hamiltonian approach. Of the four independent color schemes, 

two are investigated, the (gc)i(cg)i singlet-singlet (molecule) and the {qc)$(qc)$ 
triplet-triplet (diquark), for selected J PC states using a variational method. The 
predicted masses of triplet-triplet tetraquarks are roughly a GeV heavier than the 
singlet-singlet states. There is also an interesting flavor dependence with (qq)i(cci) 
states about half a GeV lighter than (qc)i(qc)\. The lightest 1 ++ and 1 predic- 
tions are in agreement with the observed A"(3872) and Y(4008) masses, respectively, 
suggesting they both have a ujJ ftp rather than D*D* molecular type structure. This 
is consistent with observed X, Y decays to ir^J/ip and not D*D* . 
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1 Introduction 

> 

X 

With the recent observed Higgs boson candidate at the LHC, the weak sector 
of the standard model may be approaching closure. However the strong inter- 
action component is far from complete, especially hadronic structure where 
intense searches are under way for unconventional (non qq or qqq) states. 
Indeed QCD allows for other color singlet combinations, such as glueballs, hy- 
brid mesons, tetraquarks, pentaquarks and dibaryons. As early as 1977 qqqq 
tetraquark states were proposed [1] and today there are now good candidates, 
the 7Ti(1400) and tti(1600) [2|[3], listed in the PDG meson summary table @|. 
These states have unconventional quantum numbers J PC = 1 h which are not 
possible for a qq system and have been reasonably described [3] as (qq)i(qq)i 
color singlet-singlet states referred to as molecules. In the heavy quark sector 
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no states with unconventional quantum numbers have yet to be confirmed. 
However, there are hidden charm states, some also listed in the PDG sum- 
mary, which can not be described as cc mesons. These X, Y and Z particles 
are charmonium-like but are not consistent with any cc spectrum predictions. 
The lightest is X(3872) which was discovered by Belle in 2003 [6] and has 
a narrow peak near 3872 MeV in the 7i + 7i^J/ip invariant mass distribution 
from B~ — > K~7t + tt~ J/ip decay. This discovery was confirmed by BaBar in 
the same decay process [7]. Subsequently, additional X, Y and Z resonances 
were found in variest experiments and those with confirmed J PC are listed in 
Table [U 
Table 1 



Properties of X, Y and Z mesons with established J PC . 



State 


Mass (MeV) 


T(MeV) 


J PC 


Decay 


Production B — > 


X(3872) 


3871.4 ±0.6 


< 2.3 


1++ 
or 2-+ 


TC + ir~J/vp 
jj/ip 


KX(m2) 
pp 


Z(3930) 


3929 ± 5 


29 ± 10 


2++ 


DD 


77 


Y(4008) 


4008jlg 


ZZD_ 80 


1— 


tt + it~ J/tp 


e + e~ 


Y(4260) 


4260 ± 12 


83 ±22 


1— 


ir + ir~ J '/ip 


e + e~ 


Y(4360) 


4361 ± 13 


74 ± 18 


1— 


ir + ir~i{j' 


e + e~ 


Y(4660) 


4664 ± 12 


48 ± 15 


\— 




e + e~ 



Theoretically, there have been many X, Y and Z investigations using a va- 
riety of methods, such as perturbative NRQCD 09], lattice QCD [TOlfTT] . 
effective field theory [T2ll3fl4j . QCD sum rule [T5fT6] and potential models 
[T7lll8fl9ir2T)] . all reviewed in detail elsewhere [2"Tf2"2"] . One plausible conjec- 
ture [2"3~|24f25| is these states are tetraquarks containing light, q = u or d, 
and charm, c, quarks in different color schemes. In this paper we investigate 
this further and consider two color combinations, the somewhat conventional 
molecular (qc)\(cq)\ singlet-singlet [22J and the more exotic (qc)^(qc)z triplet- 
triplet (diquark) [26] . We reserve the term exotic for the latter since it entails 
quarks in intermediate color states that are not singlets. 

Our work is an extension of a previous light tetraquark study [5] utilizing 
the QCD Coulomb gauge [CG] model, first implemented to predict a glueball 
spectrum |27,28j that was in good agreement with lattice QCD data. In this 
approach the exact QCD Hamiltonian in the Coulomb gauge is approximated 
by an effective field theoretical relativistic Hamiltonian. The bare parton (cur- 
rent quark and gluon) field operators are dressed by a Bardeen, Cooper and 
Schrieffer [BCS] rotation and variational ground state (vacuum) minimiza- 
tion. This generates a non-trivial vacuum in which chiral symmetry is dynam- 
ically broken and quark/gluon constituent masses and condensates emerge. 
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The hadrons are represented as quasiparticle excitations using many-body 
techniques such as Tamm-Dancoff and random phase approximations. This 
method was subsequently applied to mesons [2"9|3U|I31] . hybrids [3"2l33|34ll3"o] 
and light tetraquark states [5|36] . 

This paper is organized into five sections. Section [2] specifies the CG model 
which is then applied to hidden charm tetraquark states in section [3| Numer- 
ical results are presented and discussed in section [| followed by a summary, 
section |5j detailing conclusions and future work. 



2 The QCD Coulomb gauge model 



The CG model provides a comprehensive, systematic approach to hadron 
structure. It is applicable to both quarks and gluons, light and heavy me- 
son and baryon ground and excited states for any flavor and exotic systems 
involving different combinations of quarks and gluons in various color schemes. 
It also permits consistent Hamiltonian dynamical mixing between states of en- 
tirely different dressed partons while providing insight characteristic of a wave 
function picture. As further discussed below there are additional attractive 
theoretical features and it is significant to note that there are no free model 
parameters as the two dynamical constants, the string tension a and Coulomb 
interaction a s , are predetermined from the literature. 

The exact QCD Hamiltonian in the Coulomb gauge [37] is 



H, 



QCD 



Ha + H a + H, 



H, 



c 



H q = J dx^(x)[-ia ■ V + /3m]*(x) 
H g = - J(h^J- l IT(-K) -JTr(x) + B a (x) ■ B a (x) 
H qg = gj dx J a (x)- A a (x) 

J dxdyp a (x)J- 1 K ab (x,y)Jp b (y) , 
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(1) 
(2) 

(3) 
(4) 

(5) 



where g is the QCD coupling constant, \l/ is the quark field with current quark 
mass m, A a = (A a , Aq) are the gluon fields satisfying the Coulomb (transverse) 
gauge condition, V ■ A a = 0, (a = 1,2,. ..8), IT 1 = — E" r are the conjugate 
momenta and 

EL = -A a + g{l - V- 2 VV-)f abc A b A c (6) 
1 



B° 



V x A a + -gf abc A b x A c , 
2 



(7) 
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are the non-abelian chromo dynamic fields. The color densities p a (x) and quark 
currents J a (x) are 

p a (x) = ^(x)T a ^(x) + / a6c A fe (x) ■ n c (x) (8) 



J a (x) = ftfxjarifx), (9) 

where T a = ^ and f abc are the SU(3) color matrices and structure con- 
stants, respectively. The Faddeev-Popov determinant, J = det(A^), of the 
matrix Ai = V • D, with covariant derivative D afe = <5 a6 V — gf abc A c , is a 
measure of the gauge manifold curvature and the kernel in Eq. (5) is given 
by i^ a6 (x, y) = (x, a\A4~~ 1 'V 2 A4~ 1 \y, b). The Coulomb gauge Hamiltonian is 
renormalizable, permits resolution of the Gribov problem, preserves rotational 
invariance, avoids spurious retardation corrections, aids identification of dom- 
inant, low energy potentials and does not introduce unphysical degrees of 
freedom (ghosts). 

The bare parton fields have the following normal mode expansions (bare quark 
spinors u, v, helicity, A = ±1, and color vectors ^=1,2,3) 



*(*) = /(^K(k)Mk) + ^(-k)4 C (-k)]e ik ' x e c (10) 
r rlk 1 

A» = / ^ ^[aW + a«(-k)]e*" (11) 

IP(aO = -i /(^l\/|[a a (fc) " a«(-fc)]e*". (12) 



Our model's starting point is the Coulomb gauge QCD Hamiltonian. In this 
gauge, the color form of Gauss's law, which is essential for confinement, is sat- 
isfied exactly. The unphysical longitudinal degrees of freedom of gluon fields 
can be completely eliminated by solving Gauss's law. We then make the fol- 
lowing approximations: 1) replace the exact Coulomb kernel with a calculable 
confining potential; 2) use the lowest order, unit value for the Faddeev-Popov 
determinant. Therefore, He in the CG model Hamiltonian becomes 



H c G = ~\j rfxrfyp a (x)y(|x - y|)p a (y). (13) 

The confining and leading canonical interaction is represented by a Cornell 
type potential V(r) = —a s /r + ar. Previous studies with this interaction were 
in good agreement with both lattice glueball masses [27] and the observed me- 
son spectrum [30]. Performing a fourier transform, the potential in momentum 
space is 

VW) = - % (14) 
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Next, hadron states are expressed as BCS vacuum \Q) excitations involving 
dressed quark (antiquark) Fock operators B\ c , D\ C related to the bare oper- 
ators b\c, d\ c by the rotation 



5 AC (k) = cos^ ) 6 AC (k) - Asin^4 c (-k) (15) 

D AC (-k) = cos^ AC (-k) + Asin®6t c (k), (16) 

(17) 

where 9{k) is the BCS angle. Correspondingly, the dressed Dirac spinors are 



( 



1 + sin0(A;)x A 
1 — sin (p{k)(T ■ kxx 



1 / — \/l — sin (b(k)cr ■ kv A , 
V A (-k) = -= V TLjL XX \, (19) 

y ^l + sin0(A;)xA 

where 4>(k) is the gap angle from the gap equation that minimizes the energy 
of the BCS vacuum |34j and is related to 6{k) by tan((p(k) — 0{k)) = m/k. 

The four model parameters have all been predetermined [31]. The light and 
charm bare quark masses are 5 MeV and 1350 MeV, respectively, while the 
two dynamic constants are o = 0.18 GeV 2 and a s = |- = 0.4. This few pa- 
rameter, constrained model has successfully described meson, hybrid, glueball 
and light tetraquark systems yielding results agreeing with experimental and 
lattice data. 



3 Application to hidden charm tetraquark states 

For this heavier tetraquark system, the quark (antiquark) cm momenta are 
ki=i,2,3,4 and the following wave function ansatz is adopted 

\V JPC ) = J dk 1 dk 2 dk 3 dk 4 (J(k a + k 2 + k 3 + k 4 )$(f£ A3A4 (k i ) 

x^C3c:5i 1 c 1 (kiXc 2 (k 2 Xc 3 ( k 3Xc 4 (k4)|fi). (20) 

The expression for the matrix Rclcl depends on the specific color scheme se- 
lected. As depicted in Fig. [TJ the SU C (3) color algebra for four quarks produces 
81 color states, 3 <g) 3 <g) 3 <g) 3 = 27© 10® 10© 8® 8© 8© 8© 1© 1 of which two 
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Fig. 1. Four tetraquark color schemes. One is a singlet-singlet molecule while the 
other three are more exotic atoms (octet and two diquark schemes). 

are overall color singlets that can be obtained in four different schemes (only 
two are independent) depending on the intermediate coupling of two quarks: 
singlet-singlet (molecule), octet-octet and two diquark schemes triplet-triplet 
and sextet-sextet. In this article we focus on singlet- singlet (molecule) and 
triplet-triplet (diquark) schemes, which provide the lightest states. 

For the diquark configuration, the wave function is 

$£w 4 (ki) = F JPC (k l )^ 1 -X 1 X 2 \S A X A )( 1 - 1 -X 3 X 4 \S B X B )(S A S B X A X B \SX) 
(L A L B m A m B \lm) (lLjmmj\LM) (SLXM\J Mj)Y^ (k^) (k#) Y^ 7 (kj), 

where Yj^ik) is the spherical harmonic function, while S A (b), L A (b) and Li 
are the A(B) diquark total spin, orbit angular momentum and orbit angular 
momentum between the two diquarks, respectively. The radial wave function, 
F JPC (ki), is chosen to be a Gaussian, exp(— f- — Hr), where k^, kg and 

k 7 are respectively k A = k B = k 7 = - and a A , a B 

and ai are the variational parameters. 

The variational mass is then given by 



(^,Jpc\ h cg\^jpc\ 
MjPC = (yJPC\yJPC) = Msel f + M "i + M « + M ™ + Mannih ' 

where the subscripts indicate the contribution from the single particle quark 
self energy and two particle qq, qq, qq scattering and qq annihilation, respec- 
tively. All contributions to the tetraquark mass are diagramed in Fig. [2] 

Using the above specified Hamiltonian and wave function, the self energy term 
can be reduced to 
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M. 



:etra-quark — 



+ 



+ 



+ 



+ 



+ 



Fig. 2. Equal-time diagrams for tetraquarks. The first four diagrams are for the 
kinetic and self energy. The last diagram is for the annihilation term and the others 
are for the qq, qq and qq interactions. 



f dk 2 dk 3 dk 4 t+jpc /i \ 

Mself = hl (2^3(2^3(2^)3(2^3^^4 (ki , k 2 , k 3 , k 4 ) 



,., (2tt) 3 (2tt) 3 (2tt) 3 (2tt) ; 
where the single particle energy is given by 



(21) 



E i {k i )=m i sin9{k i ) + k.cosO^) - ^ J (\q\) [sin 9 (ki) sin 0(q) 

+ cos0(ki)cos0(q)q-ki]. (22) 
Similarly, the two particle contributions can be evaluated yielding 



M qq = - 



dk^ dks dkj dq 



■z^ 1 (ki)WA 1 (k 1 )z^(i4)w Aa (k 2 ) 



(2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 
x V(\q\)$ KS 3 A 4 (ki, k 2 , k 3 , k 4 )$ 3 [f A ( 7 A , A , (k' x , k' 2 , k' 3 , ki), 



(23) 



<?<2 ' 



1 /■ dk A dk B dkj rfq t t 

3 7 (27F W (2^F (kl)WAl (kl) Vls (k3) ^ (k3) 

x n| q |)$ KrS 3 A 4 (ki, k 2 , k 3 , k 4 )^ A , A , (k;, k 2 , k 3 , K), 



(24) 



2 /■ rfk A rfk B dk 7 dq t t 
-3i (27F(27F(27F(^ V V k i)^^)^ 

x V(\q\)Q Kr 2 A 3 A 4 (ki, k 2 , k 3 , k 4 )$(f A ? A , A , (k' 1; k' 2 , k' 3 , ki), 



(25) 



where q is the momentum transfer k' x — k x . Finally, the annihilation term is 
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i7 (2^(2^(2^(2^ >< *<MWWI*«Wt«> 

xv(|q||)<S,j.(k I ,fe, k 3 , k 4 )4. A 'f A ? A , A , (k;, 14, k' 3 , k' 4 ), (26) 

where now q is the annihilation momentum kx + k 4 . For the molecular case, 
each qq pair forms a subgroup and the corresponding wave function and mass 
formula are similar. The 12-dimension integrals are performed numerically 
using the Monte-Carlo method. 



4 Numerical Results 



We calculated all possible states with Si = 0,1 and/or Lj = 0, 1 and the 
numerical results for corresponding J PC are listed in Table [2j For model com- 
pleteness and further comparative insight, we also list our ccg hybrid predic- 
tions from an earlier publication [31] . It is clear that the diquark composition 
is heavier than the corresponding molecular state. Even more significant is the 
sensitivity to arrangement of flavor as (qq)\(cei) states are much lighter than 
{qc)i(qc)i for the same Li, Si configuration. The CG model therefore predicts 
that the lightest tetraquark states are molecules having quarks with the same 
flavor in a color singlet (i.e. a r\r\ c or uJ/ip type flavor structure depending 
upon spin). 

For the s wave triplet-triplet tetraquark states, the mass changes little with 
different spin. This is because there is a heavy charm quark in each color triplet 



Table 2 

Predicted spectrum for tetraquarks (singlet-singlet and triplet-triplet) and hybrid 
mesons in MeV. Unless explicitly stated, all Li and Si are zero. 



Li, Si 


J PC 


(qq)i{cc)i 




(qc)x{cq)x 






(cc)sff 


Li,Si = 


0++ 


3540 


0++ 


4129 


4428 


0++ 


3945 


Lj = l 


1-+ 


3540 


1~ 


4129 


4741 


1-+ 


4020 


L A = 1 


1~ 


4186 


1"+ 


4563 


4826 


1"+ 


4155 


L B = 1 


1~ 


4026 


1"+ 


4565 


4825 


1-+ 


4565 


5a = 1 


1+- 


3781 


1++ 


4137 


4430 


1++ 


4100 


5s = 1 


1+- 


3538 


1++ 


4137 


4430 


1+- 


3830 


Sa, Sb = 1 


0++ 


3783 


0++ 


4144 


4405 


0~ 


4020 


5a, Sb = 1 


1++ 


3783 


1++ 


4144 


4419 


3-+ 


4615 


Sa, Sb = I 


2++ 


3783 


2++ 


4144 


4444 


2++ 


3965 



s 



and the spin splitting decreases with heavier quark mass. For orbit angular 
momentum 1, the triplet-triplet tetraquark mass increases about 300 to 400 
MeV and all p wave tetraquarks have mass above 4400 MeV. 

Similarly for the (qc)i(cq)i states, the spin splitting is very small due to the 
large charm quark mass in each color singlet, and again if either La or Lb 
changes from to 1 the mass increases, now about 430 MeV . For the (qq)i{cc)\ 
flavor arrangement, the mass does not depend on the cc spin but increases 
about 240 MeV if the spin of light quark subgroup changes from to 1. The 
mass increase due to larger orbit angular momentum is also flavor sensitive, 
since for qq p waves it is about 650 MeV, while for cc p waves it is 490 MeV. 

In Fig. |H we compare predicted 1 ++ and 1 tetraquark masses to observed 
X and Y states. Only about half of the states have established J PC and to aid 
determining quantum numbers it is noteworthy that the observed F(4140), 
X(4160), F(4274) and X(4350) are comparable to the heaviest model 1++ 
states. If this identification proves correct, then the model predicts the Y (4140) 
and X(4160) are molecular states consisting of two spin clusters having a 
DD flavor composition while the V(4274) and X(4350) would be exotic spin 
1 triplet-triplet states with a diquark composition. 

Regarding the established J PC states, the predicted (qq)i(cc)i 1 ++ mass is 
3783 MeV about 3% less than observed, suggesting the X(3872) may be a 
molecular state with a ujJ /ip flavor composition. Four 1 tetraquark states 
are predicted with the lightest mass at 4026 MeV which is close to y(4008) 
discovered by Belle [38J. The masses of the next two heavier molecular states 



PDG 

Model 

r(4660) 



H4360) X(4350) 

r<4260) H4274) 

X(4160) 

7(4140) 



y<4008) 

X(3940) 

X(3915) 

X(3872) 



1— 1 ++ ? 

Fig. 3. Theoretical results (dashed lines) for 1 ++ and 1 tetraquark states com- 
pared to data (solid lines). 



4400 



> 

S 4200 



Kino 
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are 4129 and 4186, which are comparable with Y(4260) and Y(4360). The 
states with masses 4026 MeV and 4186 MeV have a 7777c flavor arrangement 
while the state at 4129 MeV is a DD type molecule. The calculated mass of 
the triplet-triplet 1 state is 4741 MeV which is also close to the F(4660) 
mass. This state would be an exotic diquark if this assignment prevails. Lastly, 
for the observed 2 ++ Z(3930) listed in Table [TJ the model tetraquark states 
are either above or below this mass by over 150 MeV, however the predicted 
ccg hybrid mass of 3965 MeV (see Table [2]) agrees quite well. Similarly the 
++ hybrid at 3945 MeV is very close to the X(3915). This state does not 
have established quantum numbers other than the G = +1 parity which is 
also the same as our hybrid prediction. The hybrid results are quite firm since 
they include both hyperfine and the non-linear (non-Abelian) component of 
the color magnetic fields (Eq. (7)). 

Finally, it is interesting that the decay channels summarized in Table 1 are 
consistent with the Coulomb gauge model flavor structure predictions for those 
observed states. Specifically, for the predicted molecular state assignments in- 
volving the X(3872), Y(4008) and Y(4260) having a tuJ/ip flavor structure, 
their decays would favor the observed tt J/ip channel. For the diquark model 
states, detailed calculations would be necessary to make firm predictions. Sim- 
ilarly, for the undetermined J PC states predicted to be DD molecules, F(4140) 
and X(4160), their decays would would favor the DD channel and, if observed, 
would be an identifying signature. 



5 Summary 

We have applied the Coulomb gauge model to tetraquark systems with hid- 
den charm having two different color compositions, singlet-singlet and triplet- 
triplet. Generally, triplet-triplet states are the heaviest for the same spin or 
orbit angular momentum. For color singlet- singlet states, we further investi- 
gated two flavor groupings corresponding to 7777,, or J/ip versus DD or D*D* 
type molecules (meson-meson states) and found that the same flavor color sin- 
glets produce a lower mass. Also, the mass from the spin part decreases rapidly 
with increasing quark mass, while the orbit angular momentum contribution 
decreases slowly with increasing quark mass. 

Our key finding is the X(3872) mass is close to the predicted 1 ++ (qq)i(cc)i 
state suggesting it is a conventional uJ /ip type molecule and not a more exotic 
diquark tetraquark state. For 1 , four states are predicted, three molecular 
and one triplet-triplet. Their masses are comparable with y(4008), F(4260), 
F(4360) and Y(4660). The masses of other states with undetermined J PC are 
close to the calculated 1 ++ results. Therefore, the J PC number of these states 
could be 1 ++ . 
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Future work will address several issues. The spin dependence will be further 
studied now including the hyperfine interaction [3T]. Also, a dynamic mixing 
analysis will be performed involving qq and qqqq states having various flavor 
distributions as well as mixing with glueball and hybrid meson exotic states. 
Finally, the CG model will be applied to b quark systems to aid experimental 
searched at higher energies. 
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